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Heavy-light and light-light weak matrix elements on the lattice * 

Claude W. Bernard 

Department of Physics, CB1005, 
Washington University, 
St. Louis, MO 63130, USA 

I review recent developments in lattice weak matrix element calculations. I focus on on /s (both with propa- 
gating quarks and in the static limit for the b quark), semi-leptonic form factors for D meson decay, form factors 
for B — > -ftr*7, and Bk- 



1. INTRODUCTION 

The main obstacle to precise determinations of 
the parameters of the Standard Model is the dif- 
ficulty in calculating hadronic matrix elements of 
weak operators. The wide range in allowed val- 
ues for hadronic quantities such as fs (the pseu- 
doscalar decay constant of the B meson), the 
form-factors for semi-leptonic B decay, and Bi^ 
(describing mixing in the K-K systems) results 
in large uncertainties in our determinations of the 
Cabibbo-Kobayashi-Maskawa angles and nit 0. 
Lattice gauge theory, at least in principle, should 
allow us to compute such matrix elements from 
first principles, i.e., with control on all sources 
of systematic error. Here, I review the current 
status of this effort. 

The talk is organized as follows. I first dis- 
cuss fg , both in the static approximation H and 
with propagating quarks. In the case of fg, the 
results show a wide spread, both among various 
groups working on the problem, and as a function 
of the lattice spacing a. I explain the reasons 
for the spread and show that the a-dependence 
seems to be understandable in perturbation the- 
ory. The variation among groups is likely due 
to the poor signal-to-noise properties of static- 
light propagators. This can introduce spurious 
dependences on the details of how the signal is ex- 
tracted. I discuss in particular the quality of the 
data obtained with different types of "smearing." 
In the case of fs with propagating quarks, the re- 
sults (at least the raw data) from different groups 



seem rather better in agreement. The main is- 
sues which I consider are the proper normaliza- 
tion of quarks with ma ^ 1, the consistency be- 
tween Wilson quarks and Sheikholeslami-Wohlert 
(SW) improved quarks, and the consistency be- 
tween the propagating and static decay constants. 

Recent results for the form factors for semi- 
leptonic decays of D mesons are then discussed 
1^. I focus on the normalization factors for the 
local and point-split currents and describe how 
their mass dependence may be understood, at 
least qualitatively. This provides a way to rec- 
oncile (or at least reduce the discrepancies be- 
tween) results from different currents with differ- 
ent heavy quark masses. 

I then describe results presented at this con- 
ference for the form factors of the rare decay 
B — !■ K*j. This is a promising new way for 
the lattice to provide a handle on the Standard 
Model. The results suggest that rather stringent 
constraints on the Standard Model could be ob- 
tained with modest improvements in the lattice 
calculations and/or the experiment. 

Finally I review some exciting new develop- 
ments in the computation of Bx on the lattice. 
All systematic errors but one appear to be under 
control, and that one is a chiral loop effect which 
should be no larger than a few percent. The re- 
sult, Bk{p = 2GeV) = 0.616 ± 0.020 ± 0.017 or 
Bk = 0.825 ± 0.027 ± 0.023 |], is an important 
contribution of lattice QCD to Standard Model 
phenomenology. 



*Review presented at Lattice '93, Dallas, Oct. 12-16, 1993 
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Results for /|;*^' from various groups 
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group 




jstat (MeV) 


FNALg 


5.7 
5.9 
6.1 
6.3 


305(15) 
269(20) 
228(21) 
237(22) 


BLS|6I 


6.3 


235(20) ±21 


APEgl 


6.0 

6.0 [clover] 


350(40) ± 30 
370(40) 


UKQCDlg 


6.0 [clover] 
6.2 [clover] 


9cq+16 + 105 
^"J"J-15-14 


P-W-C0 


5.74,6.0,6.26 
(extrap.) 


230(22) ±26 


Hashimoto 11] 


6.0 [nrqcd] 


320(30) ± 60 



2. IT' 

Table |l| gives some results for Z^''** from various 
groups. Note that the values range from a high of 
370 ± 40 MeV to a low of 228 ± 21 MeV. However, 
it is misleading to compare these final answers, 
because different groups use 1 ) different values for 
Z^'^* , the lattice renormalization constant for the 
static-light axial vector current, and 2) different 
methods to set the lattice scale, a. Since what is 
really being calculated is [f ^/rnY^^^ , the result 
for /^'^* is particularly sensitive to the value used 
for a. 

In order to make a more realistic comparison, 
1 therefore fix the scale uniformly by using com- 
putations ] [T^ of the string tension, ct, and plot 
m Fig. g the "raw" value of {f^/mY^^^a 4 vs. 
ay/a. Here, "raw" means before multiplying by 
Z^''**, except that in the case of the SW action, 
1 do correct for the difference between Z^^^^ with 
the Wilson and SW actions by using perturba- 
tion theory [ p"3|p^ ] with a boosted coupling ]|T^ 
g2 = 1.77 at /3 = 6.0 and 5^ = 1.65 at /3 = 6.2. 

Two features stand out in Fig. 0. First, there 
are real disagreements between the groups. For 
the moment, I will focus on the lower tier of re- 
sults (including the FNAL and the BLS points). I 
do this not only because it includes my own work 
but also because the FNAL group uses the 
best sources (The UKQCD results are 

somewhat higher, but only because of the rela- 



Figure 1. Raw value of {Jm^y^^^a~i vs. aa^ 
from various groups; only the relative perturba- 
tive correction between SW and Wilson quarks 
is included. Statistical errors only. Some points 
have been moved slightly horizontally for clarity. 
Key: FNAL fancy square; BLS octagon; UKQCD 
(SW action) cross; P-W-C square; APE diamond; 
APE (SW action) fancy diamond; Hashimoto 
burst. For references see Table 1. 

five perturbative correction for the SW action.) 
The second feature is the strong a dependence of 
the lower results. This is puzzling at first glance. 
If the cause were the presence of higher dimen- 
sion operators (0(a) effects), one would expect 
that that the SW {0{a) improved) points would 
show significantly less scale dependence, which 
does not appear to be the case. Furthermore, the 
results for propagating quarks seem to show con- 
siderably smaller scaling violations. (Of course, 
the size of statistical errors makes definitive state- 
ments difficult.) I believe, however, that most, if 
not all, of the a dependence can actually be ex- 
plained by the "improved perturbation theory" 
ideas of Lepage and Mackenzie. 

2.1. Tadpole improved perturbation the- 
ory and the static limit 

Lepage and Mackenzie [|l^ argue that bare lat- 
tice perturbation theory is ill-behaved because 1) 
the bare lattice coupling is a poor choice of an 
expansion parameter, and 2) lattice tadpoles are 
large. They cure the first disease by using a more 
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physical coupling constant such as gv, which is 
defined in terms of the heavy-quark potential. 
The second disease is cured by summing up the 
tadpoles with a mean field theory approach: they 
show that the simple replacement of the link U 
by its mean value mq accounts for the bulk of 
the perturbative corrections to many quantities. 
For example, for Wilson fermions, the replace- 
ment [/ — > uq just changes the hopping parame- 
ter K to K = kwq. The resulting free field theory 
has a critical hopping parameter kc = \ which 
would imply — -^j^, a relation which agrees 
reasonably well with simulations when one de- 



fines Uq 



_ /I 



trC/, 



plaq 



' 4 . In the following, I will 



define uq by Kc for Wilson quarks; that is, I take 



Mo 



1 - O.lOV 



(1) 



In computing some perturbative quantity, the 
prescription for "mean-field improvement" (or 
"tadpole improvement") is then to take out ex- 
plicit powers of uq (depending on the number 
of links which enter the definition of the ob- 
ject), replace them by the value of uq in the 
simulation (1/(8kc)), and remove the appropri- 
ate power of 1 — 0.109(7'^ -t- . . . from the remaining 
perturbative expansion. For example, the 1-loop 
renormalization factor for the local axial current 
with Wilson fermions in the chiral limit is 
Za — 2kc(1 — 0.1335^), where the ~ on Za im- 
plies that I include the standard factor 2kc relat- 
ing Wilson to continuum fermions. The tadpole- 
improved prescription would then be 

The factor of mq which is removed is really two 
factors of y^uo coming from wave function renor- 
malization on the external lines. The small coef- 
ficient in front of gy indicates that improved per- 
turbation theory is working well ||l9|. The scale 
(1/a) at which gy is evaluated is a typical one for 
non-tadpole diagrams; when the quadratically di- 
vergent tadpoles are left in, the scale is typically 
closer to 7r/a. 

An interesting feature of perturbative quan- 
tities in the static approximation, however, is 
that they are much less tadpole-dominated than 



the corresponding quantities with Wilson quarks. 
This can easily be understood by examining the 
static action: 

■ 7o> 



Ss^Y. ^-(^^)[^- - ul{x - d)h 



x-Oi 



(3) 



uoK-q] (4) 



From eq. (4) one can easily calculate, using a hop- 
ping expansion, the mean field theory static prop- 
agator. We get 



G,(x,2/)-4=^i-e-(i""o-')fe-*«)) 



(5) 



This shows that uq contributes to mass renormal- 
ization ([uuq^), but not wave function renormal- 
ization, since there is no overall factor in Gs- 
Equivalently, in momentum space one has 



Gsip) = (1 - Moe^P")- 



The uo factor looks 



like wave function renormalization at first glance. 
However, expanding about the pole, we let po = 
zlnuo +p'q and find Gs{pq) = (1 — e^^")-!^ which 
still has residue 1. Thus, unlike the Wilson case, 
there is no tadpole contribution to the wave func- 
tion renormalization of heavy quarks; one can also 
see this explicitly in the perturbative calculations 

Ell- 

In the case of the static-light axial current 
renormalization constant, one then takes out only 
a factor of ^/uq for the light quark, and divides by 
1 — Aj^g^. Using the known perturbative results 
||l^,|l4|, I then get, in the chiral limit for the light 
quark, 

^Zstat-(^)[l-5M9*)(.135-5)] (6) 



7SW 

■'A,sta.t 



[l-g^(g*)(.090 -J)] (7) 



where the superscripts "M^" and "S'VF" on Z^^stat 
indicate Wilson or SW light quarks, respectively, 
5 = .025 In (ams), and q* is the scale at which 
gy should be evaluated. The static renormaliza- 
tion constants are clearly not tadpole dominated, 
since the coefficients of gy in the above equations 
are relatively large (compare eq. (||)). 

It is not completely clear what scale q* to put 
into eqs. (§),(0). In cases where the scale af- 
ter removing tadpoles has been estimated using 
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the method of ref. [nsl, it is <l/a: for k^, one 
has q* — 1.03/a |15| ; for two non-relativistic 
QCD (NRQCD) masses, one has q* = 0.81/a and 
0.67/a (20). In the current case, the scale is esti- 
mated to be 2.36/a before tadpole removal j2l), 
but it may not be reduced as much as in the other 
cases by removing the tadpoles, since the tad- 
poles play a less important role. Indeed, since the 
time when this talk was given, the non-tadpole 
scale has been estimated by Hill and Hernandez 
||2^ ; their preliminary result is q* = 2.18/a. Note 
that (7|.(2.18/a) is quite large, ranging from 2.8 
at /? = 5.7 to 1.9 at /3 = 6.3. This makes the 
term gy (q*) (.135 — S) range from .26 to .23; with 
q* — 1/a, the corresponding numbers are .42 to 
.31. The higher order pcrturbative effects could 
therefore be quite large, perhaps 10% or even 
more, which implies considerable uncertainty in 

f stat 
J B ■ 

In Figs. ^ and |[ I put in the tadpole im- 
proved corrections to the raw static numbers, us- 
ing q* = 1/a and q* — 2.18/a respectively. Note 
that with q* = 1/a, the a dependence of the the 
lower tier of results is small. With q* = 2.18/a, 
there is still may be significant a dependence, but 
it is reduced from Fig. |l|, and the results are con- 
sistent with a constant if one drops the point at 
largest a {(3 — 5.7). One is now tempted to ex- 
trapolate to a = 0. With either scale, a linear 
extrapolation of the FNAL and BLS results gives 
Us/mf^^^cr-i 1.6 ± 0.2 at a = 0, which would 
imply /f « 190 ± 25 MeV. Extrapolation of 
the UKQCD (SW-action) points would produce 
almost the same answer. Clearly, however, there 
are large systematic errors in the result at this 
stage. I estimate that total systematic error as- 
sociated with the extrapolation, with the setting 
of the scale, with the uncertainty in the pcrtur- 
bative corrections, and with the extraction of the 
raw numbers themselves (see below) is roughly 
twice the statistical error. 

2.2. Smearing 

I now come to the issue of the large variation 
in the results from different groups. I believe this 
stems from the intrinsically poor signal-to-noise 
ratio for static-light propagators at large times 
|2|,||,ra. The problem has been dealt with by 



CO 



2,5 



2.0 - 



1.5 - 



1.0 



\r ir 



0.0 0.2 0.4 

aVo" 



Figure 2. Same as Fig. 1 but including the tad- 
pole improved perturbative corrections, with gy 
evaluated at scale q* = 1/a. 
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Figure 3. Same as Fig. 2 but at scale q* = 2.18/a. 



smearing the static quark (or both quarks) in the 
static-light sources. This can improve the situa- 
tion in two ways: 1) by producing a better over- 
lap with the ground state, thereby moving the 
plateau in to smaller times where the signal-to- 
noise is better, and 2) by increasing the statistics. 
The second point, which has not always been ap- 
preciated in the literature, will be discussed be- 
low. 

f^^^is computed from the asymptotic behavior 
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of the correlation functions 
GLs{t)^Y.(^\Aoix,t)xliO,Om (8) 

X 

" = "E(0|^s(x,t)4(0,0)|0) , (9) 

X 

Gssit) ^Y.^O\xs{x,t)xl{0,Om , (10) 

where Aq = ^17^75(7 is the local or almost local 
{i.e., point split) axial current (ft. and q are the 
static and light quark fields, respectively), and xs 
is the smeared meson interpolating operator: 

Xs{x,t) f{y)h{x + y,t)j5q{x,t) , (11) 

y 

with f{y) some smearing function. (I assume here 
that only the static quark is smeared, which is the 
most common practice.) 

I write " = " in eq. ^ because the two corre- 
lators, while equal when averaged over an infinite 
number of configurations, have different statisti- 
cal errors with a finite number of configurations. 
The smeared-source, local-sink, eq.(H), is consid- 
erably less noisy than the local-source, smeared- 
sink, eq.(|^). This was noticed with the first use 
of static-light smearing |2|]; a nice comparison 
of the results from the two appears in the recent 
UKQCD paper §]. 

The reason for the difference in statistics is easy 
to understand. Since the static quark propagates 
along a straight line in the time direction, only 
the static propagator from (0, 0) to (0, t) (coming 
only from terms with x = — y in eq. ^) con- 
tributes when the source is local. Independent of 
what the smearing function is, the result just de- 
pends on a single product of links in each config- 
uration, and hence has large fluctuations. When 
the source is smeared, many products of links en- 
ter into the result from each configuration, and 
the fluctuations are thereby greatly reduced. 

Note that the statistics issue is completely in- 
dependent of the the relative ground-state vs. 
excited-stated overlap for a given interpolating 
field. This means there may be a certain amount 
of tradeoff between good overlap and good statis- 
tics. For example, wall sources, in which both the 
static quark and the light quark are smeared over 



Figure 4. APE data for the effective mass 
from the smeared-smeared correlator for various 
smearing sizes, Ls- (3 — 6.0, SW action, 210 con- 
figs. 

the entire lattice, clearly have poor overlap, since 
the mean distance between the static and light 
quarks is much larger than the size of the ground 
state. However, they have potentially very good 
statistics since all time-like link products on the 
lattice are used. 

Figs. ^ and [s] show data for the Gss effec- 
tive mass and for Gls/Gss from APE with 
the SW action at (3 = 6.0. The heavy quark is 
smeared over a cube of side Ls] various values 
of Ls are shown. Gls is calculated with local- 
source, smeared-sink (eq. (^), which accounts 
the comparatively large errors in Fig. |^ (com- 
pared with, say, the UKQCD data [^), despite 
the large number of configurations (210). I find 
it difficult to see consistent, convincing plateaus 
in both Gss and Gls/Gss- Ls — 7 and Ls — 9 
are reasonable but Ls = 5 is somewhat doubtful 
in the former; while only Ls = 5 is convincing 
in the latter. One really needs plateaus for the 
same value of Ls in both correlators to be able 
to extract an answer. The APE group takes the 
variance between the Ls — 5 and the = 7 re- 
sults as an estimate of the systematic error here; 
however I am not convinced that this is a good 
estimator of the error. 

Fig. H shows the effective mass for vari- 
ous smearings in Gss from the PSI-Wuppertal- 
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Figure 5. APE data for the ratio Gls/Gss for 
various smearing sizes, L5. f3 ~ 6.0, SW action, 
210 configs. 

CERN group (P-W-C) I do not find the 

plateau, which they consider to start at t — 3, 
very convincing. Since their analysis determines 
the mass from this correlator and then uses it in 
determining the location of the plateau in Glsj 
I suspect there may be considerable systematic 
error involved here. I note, however, that the 
P-W-C group has varied the parameters of their 
analysis and does not find large variations in the 
results. A possible additional source of system- 
atic error arises from the rather small size of the 
lattice here (18^ x 48 at /3 = 6.26). 

Some of the difhculty in extracting /^*^* from 
simulations is illustrated in fig. ^ from our group 
^. The data, at (3 — 6.0 with cube smearings 
of side Ls = 9, is fit in two different ways: The 
dashed lines show simultaneous fits to Gss and 
Gls in the range t = (5, 10); while the solid lines 
show simultaneous fits to Gss in the range (3, 7) 
and Gls in the range (9, 13). While both sets of 
fits appear, at first glance, reasonably convincing, 
they give vastly different results. Indeed f^^'^ ex- 
tracted from the first (dashed- line) fits is « 1.5 
times larger than that extracted from the second 
(solid-line) fits. For this reason we have not felt 
confident of quoting a result for f^^^ at /3 = 6.0. 
Note that, although the number of configurations 
is small (8), the statistical errors in the data with 
our sources are at least as small as those in the 



Figure 6. P-W-C data for the effective mass from 
Gss, for various smearings. 18"^ x 48, /3 = 6.26, 
flight = .1492, 43 configs. 

data of the previously discussed groups. 

Fig. H shows our data |^ with wall sources 
a.t /3 — 6.3. Note that the errors in Gls are 
very small, but that a two-exponential fit is re- 
quired because the relative ground-state overlap 
is poor. For Gss, the errors are much larger, be- 
cause smearing at the sink does not buy anything 
in terms of statistics, but just adds noise com- 
ing from large separations between the quarks. 
These features are all consistent with the discus- 
sion above. 

Our data with cube sources at /3 = 6.3 are 
shown in Fig. |^. The errors in Gls are larger 
than with the wall sources, but the ground state 
overlap is better, and a good plateau is seen. The 
plateau for Gll is less convincing, but it does 
have a mass which is consistent with that of the 
other channel (the the channels are fit simulta- 
neously). We feel safer with the cube sources 
than with the walls because we can use single- 
exponential fits, and we extract our final results 
from cubes of side Ls = 15. Note, however, that 
the wall sources, as well as a range of different 
cube sizes, give consistent results. This data set 
is also stable under change in fitting intervals; un- 
like the f3 = 6.0 case, such shifts change the re- 
sults by only sa 3%. 

Data from the FNAL group (5| is shown in 
Fig. The sources used here are close to ideal: 
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Figure 7. BLS data for the effective mass from 
Gss (diamonds) and Gls squares. x39, f3 — 



6.0, Kiight = .155, 8 configs. The diamonds have 
been displaced by +0.4 for clarity. The dashed 
lines and solid lines are different simultaneous fits 
to the correlators; see text. 



Figure 8. BLS effective masses with wall sources 
at /3 = 6.3, X 55, 20 configs. (40 sources), 
Kiight — .150. Diamonds: Gss (displaced +0.2 
for clarity); squares: Gls- A simultaneous two- 
exponential fit to both channels is used. The 
range is (5, 13) for Gss; (3, 12) for Gls- 



their multistate fitting method (l^JT^ gives a 
source which is nearly the ground state wave- 
function. Despite this, however, it appears that 
the results could be changed appreciably by mov- 
ing the fitting intervals. Admittedly, this is not 
a "typical" but more like a "worst case" exam- 
ple. However, it should serve to emphasize that 
important systematic errors associated with ex- 
tracting the ground state exist in the data of all 
the groups. It will probably take an order of mag- 
nitude better statistics (and good sources) to be 
really confident that this systematic is controlled. 

Two new developments discussed at this confer- 
ence show some promise. First, Draper and Mc- 
Neile have devised a method of constructing 
static-light meson sources which is complemen- 
tary to the FNAL approach but produces simi- 
larly near-ideal sources. They show some prelim- 
inary data with very nice plateaus and are also 
able to extract cleanly the wave functions for ex- 
cited as well as ground states. It will be interest- 
ing to see how well they can determine /^*^* . Sec- 
ond, two groups |^,|ll[ are using NRQCD for the 
heavy quark in the heavy-light system. They seek 
to take advantage of the improved signal-to-noise 
ratio Q for NR-hght, as opposed to static- light. 



mesons. The results are quite encouraging: their 
effective mass plots are clearly better than static- 
light ones with comparable sources. Hashimoto 
evaluates of f^'^^ (see Table P; Davies et al. ex- 
tract Jb itself (see below). 

The NRQCD group ||l has also carried 
through a very nice calculation of the 6-quark 
mass. They get rrib = 4.7 ± 0.1 GeV. This has 
an important implication for the static theory. 
It has been suggested that the l/mg {Q is 
a generic heavy quark) corrections to the static 
limit would not be computable due to the pres- 
ence of power-law divergences, which could in- 
duce large non-perturbative effects. The success 
of the NRQCD calculation of m;,, which involves 
perturbative calculations of such power-law diver- 
gences and has several cross-checks, indicates that 
the concerns expressed in ref. | p9| are not likely 
to be a problem in practice. 

I have one final remark on f^^^. Given the is- 
sues relating to smearing and a dependence which 
were discussed above, I was dismayed to find the 
following comment made in ref. : 

"The result of [BLS — ref. |^] is much lower than 
all other results obtained at a fixed value of a. A 
possible explanation may be found in the fact that 
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Figure 9. Same as Fig. 8 but the sources are cubes 
(is = 13) and a simultaneous single exponential 
fit in the range (9, 16) is used. 

in [ref. Q] the results were obtained using very large 
smearings (or the waU source), which the authors be- 
Ueve to be the most coupled to the lightest state, in 
contrast with the findings of the present work and of 
[the FNAL group — ref. [|o|]." 

We have clearly stated, at least three times 
, the trade-ofF between ground-state over- 
lap and statistics which motivate our use of wall 
sources or "very large" smearings (which appear 
so large only when compared to the small smear- 
ings necessitated by the use of eq. (||) in 0). 

3. /b with propagating quarks 

Since to approach the h quark on currently 
available lattices one needs quark masses with 
ma ~ 1, it is necessary to study how lattice 
quarks behave when the graininess of the lattice is 
comparable to their Compton wavelength. Kron- 
feld and Mackenzie |32 3^ have initiated such a 
program. For zero momentum quarks the idea is 
quite simple. In a free theory the p = Wilson 
propagator is easily shown to be: 



Glatt (0 



1 



1 



-70 



1 



(12) 

where the "pole mass" M\ is related to the La- 
grangian mass mo by aMi = ln(l -|- amo), with 



"Gcont(^) , 



Figure 10. Effective masses for Gls (above) and 
Gss from the FNAL group. /3 = 5.9, 16^ x 32, 
Kiight — -157, 48 configs. 



amp = 1/(2k) — 4 for r = 1, and Gcont is the 
continuum p = propagator. The factor 2Ke"*^i 
can also be written as 1 — 6k, and it is therefore 
clear that the correct normalization for a massive 

/2k, 



free quark field is \/\ — 6k, rather than the 
which is only valid near the massless limit. 

Interactions can easily be taken into account in 
the limit of tadpole dominance. In that case, one 
just replaces k by k = kwq — k/(8kc) (see dis- 
cussion around eq. (|l|)). In the tadpole approx- 
imation, one thus has a normalization factor of 
Vl — 6k for each Wilson quark field in the limit 
K {aMi oo). 

What are the corrections? First of all there 
are ones of 0{p^). These are small if Mi :§> 
which should be the case for heavy-light mesons 
(for which \p\ ^ Aqcd ) on current lattices (for 
which aMi ~ 1 implies Mi 3> Aqcd)- A sys- 
tematic improvement procedure, similar to that 
of NRQCD, can then take these higher order cor- 
rections in to account. Second, there are 0{g^) 
corrections. These are small if tadpole dominance 
is a good approximation. Some early perturbative 

However, recall 



results seem to bear this out |33 



that tadpole dominance does not work well in the 
static limit, so one expects some fairly large per- 
turbative corrections 20% to 30%) to show up 
somewhere. 

The SW action is a systematic improvement 
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on the Wilson action, eliminating terms of 0{a) 
(but not 0{g^a)). One therefore expects it to in- 
clude the 0{aAIi) difference between the massless 
normalization {V2k) and the large mass normal- 
ization {V2Ke^). This is indeed the case. In 
the usual SW implementation, one "rotates" the 
quark field by 

+ + (14) 

One therefore automatically includes a factor 
which equals Ve°^ up to ©((amo)^). However, 
my suggestion to those using the SW action would 
be to use the exact (5 ^ 0, k — s- 0) normalization 
instead. This is "equivalent" in that its difference 
with what they are doing already is higher order 
(0(0^)), but it guarantees that the static limit 
will be recovered in the Mi — > 00 limit (up to 

Two groups have presented evidence that 
the use of the Vl — 6k norm brings the Wilson 
fermion results into closer agreement with the 
SW results, as expected. For example. Fig. |l^ 
§ shows fpVW {P is a generic heavy-light 
pseudoscalar) for Wilson quarks with the Vl — 6k 
normalization, and with the ("naive") \/2k nor- 
malization, in comparison with the SW results. 
The naive results vanish exponentially as Mp — > 
00, as expected from eq. (p^, and therefore can- 
not be consistent with the static theory. The 
Vl — 6k results are in good agreement with the 
SW results over the whole range of meson masses 
examined; this is actually somewhat surprising, 
since the two are only supposed to agree up to 
terms of O((amo)^). On the other hand, the 
P-W-C group has presented results showing 
that /p is less dependent on the lattice spacing 
if the naive, rather than the ^/l — 6k, normal- 
ization is used, although both extrapolate to the 
same value as a — s- 0. I do not understand the 
reason for these results. It is possible that the 
naive normalization may accidentally scale better 
in some limited range of Afp, but it is clear, as 
Fig. [1^ illustrates, that the naive normalization is 
completely incorrect as Mp 00, and must not 
scale with a in that limit. 

With the — 6k normalization, the results 
with propagating Wilson quarks seem to be in 



Figure 11. Comparison, from the UKQCD group, 
of results for fpMp vs. 1/M. Key: SW action 
(circle); Wilson action with (2k) i normalization 
(diamond); Wilson action with (1 — 6k) 2 normal- 
ization (square). 

good agreement with the static result. Fig. [ij 
shows the comparison between the two methods 
from our calculations at (3 = 6.3 The fit to 
both results has an excellent x^- Fig. |l^ is a sim- 
ilar comparison from UKQCD j|] at /? = 6.0 with 
SW propagating quarks. Here the solid line fit 
to the open circles and the static point looks less 
good, but, because of the correlations between 
the propagating quark points, it is actually not 
too bad: x^/dof = 1.5. The fit would improve if 
tadpole-improved perturbation theory were used 
for the rcnormalization factors Za of both the 
static and propagating axial currents. For SW 
propagating quarks, Za actually is larger after 
tadpole improvement; while in the static case it 
is reduced. In Fig. |l^, I have attempted to show 
how the points would move after tadpole improve- 
ment. There is some uncertainty here since it is 
not completely clear how to apply tadpole im- 
provement in the SW case {e.g., whether to use 
the Wilson or SW Kc — there is considerable dif- 
ference). However it is clear that the agreement 
between static and propagating will get better. 
Use of tadpole improvement should also reduce 
the large uncertainty in scale in the UKQCD re- 
sults, which is due to a low value of /jr. Tad- 
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Figure 12. The combined (propagating and 

static) analysis for /pM|/(l + ^ln(aMp)) at 
P — 6.3 from the BLS group. The propagating 
quarks are normahzed with (1 — 6k;) 2. The sohd 
hne is a covariant fit, quadratic in 1/Mp, to the 
marked points, and has x^/doi — 1.9/3. If the 
renormahzation were done exactly as advocated 
in sec. 2.1, the static point would move down 3% 
(using q* — 2.18/a); the propagating point, 1.5% 
(using q* = l/a). 

pole improvement will bring up the corresponding 
renormahzation constant, making it closer to the 
non-perturbative value of ^ 1.09 Note also 
that a determination of fs using non-relativistic 
heavy quarks (on the same UKQCD configura- 
tions) gives a result consistent with the other two 
methods |2^. 

In early work |36| , fits of static and prop- 
agating results seemed to show agreement be- 
tween the two approaches, even though the naive 
^/2k normalization was used in the propagat- 
ing case. However, such fits were non-covariant; 
they did not take into account the correlations in 
the propagating-quark data. In ref. |^ , we show 
that a reasonable-looking non-covariant fit can be 
made between the static point and moderately 
heavy, naively normalized, a propagating results. 
However, once correlations are included, the fit 
looks very bad and has x^/dof = 16/2. Clearly: 
(a) results with naively normalized propagating 
quarks are inconsistent with the static approach, 
and (b) covariant fits with meaningful s-re cru- 



1 

Figure 13. The combined analysis for fpMp + 
^ In(aMp)) atp = 6.0 from UKQCD. SW prop- 
agating quarks are used. The dashed line is a fit 
to the open circles and the open square is its in- 
tercept at 1/Mp = 0. The solid line is the fit 
with the static point (cross) included. The solid 
circles show my estimate for the central values if 
tadpole improvement were used. 

cial when testing the consistency of various ap- 
proaches. 

I also call attention to the recent determina- 
tions of heavy-light decays constants in full (un- 
quenched) QCD [Q. The results are generally 
consistent with the quenched results, although 
at heavy mass (close to the B) there tends to 
be some disagreement (the full QCD results are 
higher). This is not likely to be the effect of 
quenching, especially since the dynamical quark 
masses used are not very light. Rather, current 
computer limitations force the full QCD calcula- 
tions to be done at rather strong coupling. Near 
the B, ma can then be quite large, and the 0{p'^) 
corrections mentioned above can become quite 
important. In particular, the coefficient of in 
the kinetic energy is no longer 1/ (2Mi), where 
Ml is the "pole mass," but 1/(2M2) Since 
M2 ^ Ml for large Mia, the full QCD points ac- 
tually correspond to considerably heavier meson 
mass. If this effect were put in, the disagreement 
between quenched and full results would likely go 
away. Note that the correction Mi — > M2 has 
already been applied to the data in ref. 
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Table 2 

Results for fs and fo from various groups 



group 


fB (MeV) 


Jd (MeV) 


BLS|6| 


187(10) ±37 


208(9) ± 37 


UKQCD|[ 


1 






P-W-C[34 


1 


204(50) 


212(40) 


APE|^ 




230(30) 


NQRCD[27l 


160(10) 




HEMCGCI 


371 


200(10) ±48 


250(5) ±45 



Table 3 

Comparison of other BLS and UKQCD results. 
Decay constants are in MeV. 





BLS§ 


UKQCD [| 




207(9) ± 40 


194+^+;^^ 


/d. 


230(7) ±35 


212+1+^" 


/b, / fs 


1.11(2) ± .05 


1 99+. U4 
-•^•^^-.03 


Id, I fo 


1.11(2) ± .05 


i.i8t:'^^ 



Table ^ shows results from various groups for 
/b and fo, and Table ^ compares results from 
BLS and UKQCD for other decay constants and 
ratios. The results are generally in quite good 
agreement. At the current level of errors, differ- 
ences in procedure do not show up significantly 
here. For example, in the case of the BLS vs. P- 
W-C results, different normalizations of the prop- 
agating quarks {-s/l — 6k vs. V2k) tend to cancel 
the effects of different fg'^^ (see sec. 2.2), produc- 
ing quite similar final results for Jb. The only sig- 
nificant difference between the BLS and UKQCD 
results occurs in the ratios, especially fsj fa- 
This difference can be traced mainly to the fact 
that UKQCD do not include the static values for 
the ratios in their analysis: both groups get rel- 
atively low values (1.11, 1.14) in the static limit. 
Because of the problems with the static signal, it 
is not clear at this point which procedure is more 
reliable. 

A result not shown in the tables is the HEM- 
CGC value using full QCD (/3 = 5.3, Wilson 
quarks): fo, 345(5) ± 48 MeV. I do not un- 
derstand the reason for this large value compared 
to the quenched results. HEMCGC suggest that 
it may be due to the larger lattice spacing in these 
full QCD calculations. 



4. SEMILEPTONIC FORM FACTORS 

Four groups have recently obtained new results 
for D ^ K and/or D ^ K* ||-|l|. On the 
whole, the results are in agreement with previous 
work 1 ^ , ^ . Nice plots of the form factors as 
functions of (the square of the 4-momentum 
transfer) have been obtained. Unfortunately, the 
statistical errors are not in general appreciably 
smaller than in the earlier work. This is pre- 
sumably due to the inherent noisiness associated 
with lattice propagators of particles with non-zero 
three-momentum. Just as for static-light mesons 
||2^,|2^,|l^,Q , the signal-to- noise ratio for mesons 
with non-zero momentum falls exponentially be- 
cause the average of the squared propagator has 
non-zero overlap with a state which has less than 
twice the meson energy. In this case, the overlap 
is with a state of two mesons at rest. 

One interesting issue concerns the value of the 
form factor ^2(0) (or the ratio A2(0)/Ai(0)). Be- 
cause of the size of the errors, all the new and old 
results agree within 2a. The new work has, how- 
ever, has helped to change the qualitative picture. 
Initially the ELC group [|2) found an A2 (0) con- 
sistent with zero (0.19 ± 0.21 or, later, 0.4 ± 0.4) 
and with the initial experiment (0.0 ± 0.2 ± 0.1) 
On the other hand, the BES [|| result 
was distinctly different from zero (A2(0)/Ai(0) = 
0.70± 0.16to i§). The recent lattice calculations 
are now giving an ^2(0) bounded away from zero: 
A2(0)Mi(0) = 0.7±0.4 H; ^2(0) = 0.67±0.44 
or 0.72 ±0.50 (on two sets of configurations). 
A non-zero ^2(0) is at present also found by ex- 
periment: A2(0)Mi(0) = 0.82 ±0.23 ±0.11 11; 
^2(0) = 0.44 ±0.09 ||]. 

Although rcf. |3^ takes a look at the extrap- 
olation to the B meson, the other lattice studies 
to date on heavy-light semileptonic form factors 
have focussed exclusively on D mesons. The is- 
sue of normalization of fields and currents when 
amq ~ 1 is thus not as crucial here as for Jb- 
Still, it is significant numerically. 

Consider the normalization of the lattice vec- 
tor current for large am. For the local cur- 
rent, "t/|°cai^^^ _ ip{x)'^fj_ip{x), the arguments of 
refs. 1 32 1^ for normalizing the quark field would 
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just give a normalization factor 

0{g^) = 2Ke^''-'i" 



^local = 1 _ 6k 



Oig') , (15) 



where the O(g^) terms are expected to be small 
since the tadpoles are summed by using k and 
Ml. 

The normalization of the conserved current, 



T rcons 



Table 4 

Values of "Zj?'^'''" from Table 7 of ref. [43], and 
the values as calculated perturbatively using the 
ideas of refs. [31,14]. J is the heavy-light and V 
is the light-light vector current, respectively. 



matrix clem. 


ref. 11 


pert.th. 


{ViVi) 


0.659(3) 


0.70 


{KJoD) 


0.87(-) 


0.83 


{KJiD) 


0.75(1) 


0.70 


mjiD) 


0.66(2) 


0.70 



is more subtle. In the mean field approximation, 
it is not however hard to see how the current must 
be normalized in the limit am — > oo. Consider the 
matrix element 



^(oiv(f,i)v^;.(o, 0^(0, 0)10) 



(16) 



V;™'' to 



V. 



local 



and compare the case - • ^ , ^ - . ^ 

Because the quark just propagates in the time di- 
rection for large ma, the difference between V^^""'' 
to Vjj"'^'^' will simply be that the former has an 
extra factor of uo (the mean field value of the 
explicit link in its definition) but is missing a fac- 
tor of e^^^i"^ (since the point splitting allows it 
to skip one hop in the time direction). We thus 
expect 



-Mia 



Uo 



^T-local 



2k 



(17) 



Eq. (^^ is actually exact; the conserved current 
obeys a Ward identity which can be used to show 
that J2x 2kVq™^^{x, t) just counts the total charge 
@- 

The spatial components V^°'^^ are however nor- 
malized differently. Comparing eq. (|l^) for con- 
served and local currents, we see that the former 
again has an extra factor of mq but is not miss- 
ing the factor e"**^^", since no hops in the time 
direction are saved. We thus expect 



;7Cons, 



Uo 



local 
V 



— 2Ke 



Aha 



+ 0{g^) . (18) 



Note that although the /i = and — j com- 
ponents of V^°"'' enter the same Ward identity, 
this does not guarantee that they are normalized 



the same. The identity involves the divergence 
of V^°"'', rather than V'^Q'is itself, and one would 
need to take into account the difference between 
gma _ ^ g^j^j jTQj. i-j-^g discrete time derivative. 

In the literature, the quantity 



{p\vr"\^) 

(p|l/local|0) 



'7 local 



(19) 



where p is a generic vector state, is some- 
times used as a "non-perturbative definition" of 
^■locai^ However, the above arguments show that 

that the " 



is not fixed by the Ward identity and 
non-perturbative Zi?'^'*'" should in fact 



be roughly independent of the quark mass, un- 
like the correct Zy'^^'^ of eq. disl). Indeed, this 
rough mass independence is found in simulations 
1^,^,^. The independence is not evidence for 
a failure of the Kronfeld-Mackenzie ideas, as is 
suggested in refs. |^,^. On the contrary, it is 
evidence in favor of ref. 



|32| 



One can go further. Table 7 of ref. 38 gives the 
values of the "non-perturbative Zy"^** from vari- 
ous comparisons with the conserved current. Us- 
and a tadpole-improved pertur- 



ing eqs. (|l7|),(|li 

bative calculation of Zy'^^'^ with coupling gv{y/a), 
I am able to (roughly) reproduce their numbers, 
as shown in Table ^. I am neglecting the non-zero 
3-momentum that appears in the last two matrix 
elements. 

Some groups {e.g., refs. j4|,0) have compared 
with experiment their results for the vector meson 
decay constant, l//y, as a function of Mp/My 
(the square of the pseudoscalar to vector meson 
mass ratio). Such comparisons can be very mis- 
leading since it is known that Wilson quarks get 
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much too small a value for the My — Alp splitting 
for heavy quarks. Thus this is not a good way to 
test various large-am normalizations of the vector 
current. 

Finally, the normalization factor V e^^^°-, which 
is not included in the analysis of ref. , is as as 
large as 1.3 for their heaviest quark masses. Thus 
the extrapolation of form factors to the B mesons 
in that paper should be taken as qualitative only. 

5. FORM FACTORS FOR B K*j 

The penguin decay b sj [^9| provides an ex- 
cellent test of the Standard Model (SM). It is a 
short-distance decay at the quark level (the loop 
with a virtual t-quark dominates), so the inclu- 
sive rate can be calculated perturbatively. In the 
SM, the branching ratio BR{h sj) is a slowly 
varying function of the top quark mass and is es- 
sentially independent of the CKM angles in the 
3-generation case. These features make it very 
sensitive to physics beyond the SM. 

The branching ratio for the exclusive decay 
B A'*7 has been measured experimentally to 
be (4.5 ± 1.5 ± 0.9) x 10"^ @. For theory to 
make contact with this result, the ratio 



Rk* 



T{B K*j) 



(20) 



r(6 S7) ' 

is required. Phenomenological evaluations of Rk* 
vary from ~ 1% to ~ 97%, so the lattice can help 
out greatly here. 

In leading logarithmic order one only one needs 
to compute the form factors Ti{q'^), defined by 

{K*{k)\-sa^''qJ-^-^^b\B{p)) = 



where the operator on the LHS arises from the 
short distance expansion of the penguin loop, 
q = p — k is the photon momentum, and the 
definitions of the coefRcients can be found for 
example in refs. ||5l|,|5^. The form factors obey 
Ti(0) — 12(0), and T3 does not contribute for a 
real photon, so Rk' is determined by Ti{0). 

After a first look on the lattice |53|, two eval- 
uations of the form factors have been attempted 



recently ||51|,^ . An extrapolation to the B mass 
is required as usual. The two groups differ in 
how they perform the extrapolations. BHS (si] 
check that pole dominance relates Ti(0) = 12 (0) 
to r2(q^,jjx), then extrapolate T2(g^ax) to the B 
using the relation derived from heavy quark ef- 
fective theory: 

y/M^T2{ql,J^A+^ , (21) 

and finally use pole dominance to get 12(0) from 
T2(gmax) at the B. UKQCD || calculate Ti(0) 
directly on the lattice (with some mild interpola- 
tion) and then extrapolate Ti(0) linearly in 1/Mp 
to the B. 

Each method has some disadvantages: BHS use 
pole dominance over a wider range (QmaxZ-^^s ~ 
0.65) than can be checked with the current lat- 
tice data {ql^jMj,<0.5). UKQCD do a "blind" 
(i.e., without theoretical guidance) extrapolation 
of T'i(O) to the B mass. Obviously I prefer the 
BHS approach, but both methods are plausi- 
ble. The difference in the final results gives some 
measure of the systematics: BHS get ri(0) 
0.10 ± 0.01 ± 0.03; UKQCD get Ti(0) 
(assuming approximate independence of the spec- 
tator quark mass). 

Using their result for Ti(0), BHS find Rk' = 
(6.0± 1.2±3.4)%. Combined with the experimen- 
tal BR for B K*'y, this produces a la bound 
which is, unfortunately, uninteresting: mt>100 
GeV. However, because BRib 57) is a slowly 
varying function of riit, it would not take much 
improvement in the lattice and/or experimental 
results to produce a rather stringent bound. 

6. RECENT DEVELOPMENTS ON Bk 

In the past year, there has been considerable 
progress in the computation of Bk, the if"- 
K'^ mixing parameter. First of all, the cal- 
culations of the 0{g'^) perturbative corrections 
to all relevant lattice operators have been com- 
pleted 1^-^. These include 4-quark opera- 
tors of the following types: gauge-noninvariant 
(Landau gauge) "unsmeared" (2^) [||,|7|j5|] and 
"smeared" (4^)|5^,Q, and gauge- invariant (2"* 
with links) p^ , p7 |. Without the perturbative cor- 
rections, results for Bk from different operators 
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Figure 14. Comparison [60] of Bk before (open 
symbols) and after (solid symbols) renormaliza- 
tion correction. The superscript "inv" denotes 
gauge-invariant operators; "non-inv," unsmeared 
Landau-gauge operators. 



differ by as much as ~ 15% (^ 7cr) at /3 — 6.0 

f, |60|] and as much as ^ 7% (~ 5a) at P ~ 6.2 
With the perturbative corrections included, 
results with different operators agree within er- 
rors An example is shown in Fig. |lj, 
taken from ref. |60|. Note that with the bare 
lattice coupling the disagreement between differ- 
ent operators is reduced but not removed [59| ; 
full agreement occurs when one uses an improved 
("boosted") coupling a la Lepage and Macken- 
zie [|5|. (Tadpole improvement is not used here 
because various operators involve different num- 
bers of links.) The results 62 6^ also agree 
between groups. 



The second advance involves the comparison 
of quenched and full QCD results for Bk- At 
the level of the statistics (~ 2%), no difference is 
found 63 1 . The two dynamical quark masses 
in the full QCD simulations are degenerate and 
roughly equal to the average quark mass in the 
kaon (niq « TOs/2). One can easily estimate, us- 
ing the known chiral logarithms ||6^ , |65|] , how large 
a difference one would expect in the full theory 
between the simulated case and the physical sit- 
uation with nid ~ rriu ~ 0. The result is a differ- 
ence of only 3% to 4% . This of course needs 
to be checked in simulations. Unfortunately, the 



quenched theory cannot be used to estimate the 
difference between a Bk with degenerate quarks 
and one with a light rrid . Quenched chiral pertur- 
bation theory shows that the chiral logs for 
a Bk with non-degenerate quarks (in contrast to 
the degenerate case) are different in the quenched 
and full theories. 

The third advance is the understanding of the 
lattice spacing errors. At Lattice '91, the Bk data 
||6^ showed a strong dependence on a, and it was 
unclear how to extrapolate to the continuum: A 
linear extrapolation in a gave BkO-s = Bk = 
0.66(6); a quadratic, Bk = 0.78(3). This dif- 
ference dominated the systematic errors. Sharpe 
Q has now shown that the lattice spacing er- 
rors are 0{a^). The basic idea is simple: the 
staggered fermion action has corrections only at 
O(a^). So corrections of 0{a) to Bk could come 
only from the dimension-6 weak operators them- 
selves, through mixing with dimension-7 opera- 
tors. However, a detailed enumeration shows that 
no such operators exist with all the right sym- 
metries: hypercubic group, flavor symmetry, and 
individual axial rotations of each flavor. 

Once one knows the lattice spacing errors are 
0{a'^), it is straightforward to include the pertur- 
bative corrections and extrapolate existing data 
of Gupta, Kilcup and Sharpe to the contin- 
uum. The result is Q 

Ba'(NDR, 2GeV) = 0.616 ± 0.020 ± 0.017 (22) 
= 0.825 ± 0.027 ± 0.023 , (23) 

where the continuum with A— = 300McV is 

MS 

used to obtain the second equation. The only er- 
rors not included above are those due to quench- 
ing and to the use of degenerate quarks. As dis- 
cussed above, the quenched and full theories agree 
for degenerate quarks at a scale roughly that of 
the quenched /? — 6.0 (a^^ ~ 2GeV). One needs 
however to check that the theories continue to 
agree at smaller lattice spacings. In addition, 
the effects of non-degenerate quarks in the full 
theory need to be studied. However, given the 
size of the chiral logs, I would be surprised if this 
amounted to more than a 5% or 6% error. Thus 
if one at most doubles the systematic errors in 
eqs. (p^),(23), one has a result which is likely to 
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be very reliable. This is a major accomplishment 
of lattice QCD. 

I'd like to briefly mention one other result 
which has implications for light-light weak matrix 
elements. Kuramashi and collaborators [|66| have 
obtained results for tt-tt scattering in QCD. Their 
method makes possible the first calculation of the 
/ = scattering length. The same method may 
ultimately allow the computation of the AI =1/2 
K — > TTTT weak amplitude, a computation which 
has bedeviled some of us in lattice gauge theory 
over the past ten years. 

7. CONCLUSION 

The basic results have already been summa- 
rized in the introduction. I would just like to add 
the observation that we owe a great debt to Lep- 
age and Mackenzie for showing how to understand 
and control lattice perturbation theory ||l^ . This 
has made possible many of the advances discussed 
above. 
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